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Abstract 

We provide two methods of producing the Q-operator of XXZ spin chain of higher spin, one 
for iVth root-of-unity q with odd N and another for a general q, as the generahzation of those 
known in the six-vertex model. In the root-of-unity case, we discuss the functional relations 
involving the constructed Q-operator for the symmetry study of the theory. The Q-operator 
of XXZ chain of higher spin for a generic q is constructed by extending Baxter's argument in 
spin-i case for the six- vertex Q-operator. 
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1 Introduction 



The objective of this paper is to construct the (5-operator(s) of XXZ spin chain of higher spin. 
The Q-matrices were invented by Baxter to solve the eigenvalue and eigenvector problems of the 
eight-vertex model [2, 3, 4, 5]. The first one, denoted by Q72, was introduced by Baxter in 1972 [3] 
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on some special " root-of-unity" crossing parameter rj, where he computed eigenvalues of the eight- 
vertex model by establishing a functional equation, called the TQ-relation, between the transfer 
T-matrix and auxiliary Q72-matrix. One year later in the computation of eigenvectors of the transfer 
matrix, he used a different Q-ma'trix in [4]. Further later, another Q-operator was constructed in 
[5] for the study of eight-vertex model at the general crossing parameter rj. All those Q-opcrators 
satisfy the same TQ-relation, which implies the Bethe equation of the model. In recent years, the 
degeneracy of the spin-^ XXZ Hamiltonian with the extra s/2-loop-algcbra symmetry was found 
and extensively analysed in [11, 14, 15, 16, 17] when the (anisotropic) parameter q is a root of 
unity. The root-of-unity symmetry about the degenerate eigenspace of the six-vertex model has 
been further extended to the eight-vertex model [12, 13, 18] and the XXZ spin chain of higher spin 
[25, 30]. Furthermore, a set of (conjectural) functional relations in the root-of-unity eight- vertex 
model was proposed by Fabricius and McCoy in [19], as an analogy to functional relations known 
in the AT-state chiral Potts model (CPM) [9, 10], where the chiral Potts transfer matrix in theory of 
CPM corresponds to some proper Q-operator of eight-vertex model that encodes the root-of-unity 
symmetry properties. Subsequently the functional-relation aspect on the symmetry of solvable 
models suggested by Fabricius-McCoy comparison led to the discovery of the Onsagcr-algcbra 
symmetry in the superintegrable CPM, parallel to the s/2-loop-algebra symmetry in six-vertex 
model [26, 27, 28]. On the other hand, one may regard the CPM as a "model" theory in studying 
the degeneracy of solvable lattice models due to the better understanding of Baxter's Q-operator in 
CPM (i.e. the chiral Potts transfer matrix) in the functional-relation setting [1, 6, 7, 24]. Along this 
line, the investigation was successfully carried out in the root-of-unity six- and eight-vertex models 
[21, 29, 31]. However in the case of XXZ spin chain of higher spin at roots of unity, the functional 
relation study has not been conducted so far albeit the sZ2-loop-algebra symmetry of the system was 
already shown in [30] using the algebraic Bethe-ansatz method and representation theory. Indeed 
the Q-operator of XXZ spin chain of higher spin, except in the spin-i case, has not been explored 
in literature to the best of the author's knowledge. In the present article, we provide two methods 
of producing Q-operators of XXZ spin chain of higher spin, one for a root-of-unity q and another 
for the generic q, where in the root-of-unity case, the order of q, denoted by N, will assume to be 
odd. For the root-of-unity q, the Q-operator is constructed by the same mechanism of producing 
the Q-operator of the root-of-unity six- vertex model in [29] , a method imitating the construction of 
Baxter's (572-operator in [3]. For the general q, we extend the argument of producing the six- vertex 
Q-matrix in [5] to the general case of XXZ chain with higher spin, then obtain the corresponding 
Q-operator. 

This paper is organized as follows. In section 2, we briefly review some basic facts about the 
XXZ spin chain of higher spin, that are needed for later discussions. In section 3, we construct 
the Q-operator of the spin-^^ XXZ chain at Ath root-of-unity q with d < N, for the functional- 
relation study incorporating with the root-of-unity symmetry found in [30]. Section 4 contains 
another Q-operator construction of the XXZ spin chain of higher spin, but valid for a generic q. 
By techniques in [5] of producing the Q-operator of six-vertex model, we obtain a Q-operator of 
XXZ chain of higher spin with an arbitrary value of q. We close in section 5 with some concluding 



2 



remarks. 



2 The XXZ Spin Chain of Higher Spin 

We first introduce here some basic concepts about the XXZ spin chain of higher spin needed for 
later discussions. The summary will be sketchy, also serves to establish notations, (for more details, 
see e.g. [22, 30, 32] and references therein). 

In this paper, the XXZ chain of spin for positive integers d>2 will always assume with the 
even chain-size L. We use e*^ (0 < A; < d — 1) to denote the standard basis of C'', and (0 < fc < 
d— 1) the dual basis. The local spir 
and is the total-spin operator: 



d— 1) the dual basis. The local spin-operator of will be denoted by := dia[^^, . . . , — 



e=i 

The L-operator of XXZ chain of spin-^^ is the matrix with the C^-auxiliary and C'-quantum 
space: 



V Li,o Li,i J 



0<k,k'<d-l 



(2.1) 



for s G C with Lj jj^' being zeros except 

-O.Ofc 



I k I d—k—1 

Ln.n7, — Li 1 



d-k-2 



d-k-1 
I k+1 



a{q'^s) 

„k+l _ 



{0 <k <d-l), 
-''-^ {0<k<d-2). 



(2.2) 



Hereafter a(s) will denote the (g-dependent) s-function 



a{s) = sq "2 — s . 

It is well-known that the C^-operators Ljj give rise to the d-dimensional irreducible representation 
oiUgish) (see, e.g. [22, 23]), 

d-2 - — -(d-2) - — -1 

lo,o = sq 2 K 2 - s '^q 2 ^2, Lo,i = (g - g ^)e_. 



d-2 - — -1 



-(d-2) ■ 



Li,o = (g-g )e+, li,i = sq 2 K2 - s ^q 2 K2, 



with = q^^and Ke±K-^ = q^'^e±, [e+,e_] = ^gZf-i^ ■ The L-operator (2.1) satisfies the YB 
relation 

it;(^)(L(.) 1)(1 L(.')) = (1 (g) L(.'))(L(.) 1)R{^) (2.3) 



with the i?-matrix 



Ris) 



( s ^q — sq ^ 






s ^ — s q — q ^ 
q — q~^ s"^ — s 







s~^q — sq~^ j 
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The monodromy matrix of size L, 



^l,{s)=(^^'^^ ^l^M, L,(t) := L(t) at site ^, (2.4) 



e=i 



again satisfies the YB relation (2.3), hence its trace, cahed the transfer matrix of the XXZ chain 
of spin-^^, forms a commutative family of ((g) C'^)-operators: 

tis) = A{s) + D{s). (2.5) 

Then [t{s),q^''] = where q^'' =(8) . As the i?-matrix at s = g is of rank 1, the quantum 
determinant of (2.4) is defined by 

R{q){^Leiqs) 0aux 0aux ^U{s)) = 
(1 ®U{s)){^U{qs) 0aux ^)R{q) =■ detg ® U{s) ■ R{q), 



or equivalently, 



A{qs)C{s) = C{qs)A{s),B{qs)D{s) = D{qs)B{s), 

B{s)A{qs) = A{s)B{qs),C{s)D{qs) = D{s)C{qs); 
detg ® leis) = A{qs)D{s) - C{qs)B{s) = D{qs)A{s) - B{qs)C{s) 
= D{s)A{qs) - C{s)B{qs) = A{s)D{qs) - B{s)C{qs). 



(2.6) 



For the local L-operator (2.1), detgL(s) = a(s)a(g'^s), so the quantum determinant of (2.4) is equal 
to a{s)^a{q'^s)^. 



3 The (^-operator of XXZ Spin Chain of Higher Spin at Roots of 
Unity q 

In this section we study the Q-op crater of the XXZ chain of spin at the A'^th root of unity q 
with odd N for 2 < d < N . We will take q and q2 to be primitive A'^th roots of unity. 

As in the case of root-of-unity six- vertex model in [29], the Q^j-matrix of the spin-^^ XXZ 
chain is constructed from an S-operator, which is a matrix with C^-auxiliary and C^-quantum 
space, S = (Sjj)jjg2Y ~ Sjj(s)), where the C^-auxiliary space has the basis indexed by 

Z7v(:= Z/AZ), and Sij are C^-operators. The Qij-operator is defined by 

L 

Qfi = tr^jv((^S£), = S at site i. 
1=1 

Then X.Qr = tr^2 qn{<^^^i U^), where = U at the site I, and the local-operator U is the 
matrix with the ® C^-auxiliary and C'-quantum space: 

" I Li nS Li,iS / ■ 
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The operator iQr will decompose into the sum of two matrices if we can find a 2N by 2N scalar 
matrix M (independent of s) of the form 

^^{^S In)' " ' ' ^^'^^ 

such that M~^UM=( j. One can express M~^UM by 



«jgZjv 



where A(5j), C{di,6j), D{5i) are the C^-operators 

Ai5j){s) := Lo,o + \-o,iSj, D{Si){s) := -^^Lo,! + Li,i, 
C{di,6j){s) := -SiLofl + Li^o - <^iLo,i<^j + Li,!^^. 

The above operators satisfies the following commutative relations for arbitrary di,6j 



(3.2) 



(3.3) 



C{6i,6,){s)D{6i){qs) = D{6i){s)C{Si,Sj){qs). 
In fact, by the quantum determinant property (2.6) of the L-operator, one finds 

(Li,o + li,iSj){qs)A{Sj)is) = A{Sj){qs)ili,o + li,i5j){s), 
{-SiLo,o + li,o){s)DiSi){qs) = D{Si){s){-diLo,o + Li,o)(gs). 

Then the equality, C{6i, 5j) = —diA{dj) + Li^o + Li,i5j = — (5jLo,o + Li,o + D(5j)(5j, in turn yields the 
relation (3.3). 

Lemma 3.1 Let A{(3),D{a) and C(q;,/?) be the operators defined in (3.2) for a,/? G C* and q 
an arbitrary number (not necessary a root of unity). Then the criterion of a, [3 with the zero 
determinant function, detC(a, /5)(s) = 0, is: (3 = q'^~^^'^^a for some < k < d — 1. In this 
situation, C{a,f3){s) has one- dimensional kernel and cokernel for all s G C, generated by a non- 
zero vector v{= v(s)) € C^, v(= v(s)) G C^* (unique up to non-zero multiples of scalar s -functions) 
with C{a, /3){s)v = vC(a, /9)(s) = 0. Furthermore the following relations hold: 

A{P){sHs) = a{s)^y{qs), D(a)(.)v(.) = a(g'^-is)^v(|), 
v(5)A(/?)(.) = a(5'^-i.)|^v(|), v(.)D(a)(.) = a{s)^^^qs), 

v{-) 1 v{qs) 

wherev^s)^ = eo(v(s)), \/{s)'^~^ = v(s)(e'^~^), thee^- and e^-i- component ofv{s), v{s), respectively. 

Proof. The relation /3 = q^~^~'^''a with < A; < d — 1 follows from the equality 

d~i 

detC(a,/?)(s) = det(/3 - aq'^^^) [] a{q''s). 

k=0 
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By (2.2), the C^-operator C{a,/3){s) in (3.2) has the rank at least d — 1 in general. When (3 = 
qd-i-2k^^ it has the one-dimensional kernel and cokernel for all s G C, hence generated by the non- 
zero vector v(s) G C'^,v(s) G C^*, respectively. In this situation, by (3.3), C (a. /?)(g,s)(A(/?)(s)v(s)) = 
C{a,l3){q^^s){D{a){s)v{s)) = 0. Since v(gs), v(g-^s) can be characterized as the solution, unique 
up to scalar s-functions, for the equation 

C{a,P){qsHqs) = C{a, ^){q-h)y{q-h) = 0, 

the expressions of A(/3), D{a) in turn yield the first two relations in (3.7) about v. Similarly, follow 
the other two relations in (3.7) about v. □ 

Remark. The closed form of the vector v and covcctor v in the above lemma is in general hard to 
obtain except the cases P = q'^^^~^^a, where v = Ylk=o ^fe^^' ^ = Sfe=o ^k^'' given by-*^ 

^=a-fc(sg^)±^ (3.5) 

Here = is the g-binomial for integers < n < m with [n] = ^£^t-, [n]\ = Yli=i[i] 

and [0]! := 1. Indeed in the above cases, one can directly verify the relation (3.4) using the formula 
(3.5). 

In the procedure of constructing Qij-operator in the A?^th root of unity q, the necessary condition 
for Si,Sj with the zero lower blocktriangular matrix of M~-'^UM, C{Si,Sj){s)Sij = for Sjj ^ 0, 
is the vanishing determinant of C{di,6j){s) for all s. By Lemma 3.1, Sj = q'^~^'''^^5i for some 
< A; < cZ — 1, and the non-zero Sjj with C(5j, (5j)(s)Sjj(s) = is expressed by 

Si,j{= S„(s)) = Vij(s)rij, Vij(s) € C^, Tij G C''* (3.6) 

where yi,j{= Vij(s)) is the kernel vector of C{Si, Sj){s). Hereafter we shall always choose the above 
dual- vector Tij as parameters independent of s. By (3.4), we find 

D(5,)(^)S,,(.)=a(g'^-i.)^g|i;^S,,(f), 

where Vi.jis)^ = eo(vj.j(s)). For the QL-operator, we repeat the above working, replacing S, L(s)S 
by S, SL(s) with the same M. We form the product Ql^ using Ql = tr ^jv(0^=i S^), and S = 
ihj)i,jeZ^ with 

Si,j{= Sijis)) = rijvijis), %,j e v,j(s) G C''*, (3.8) 

where \/jj(= Vjj(s)) is the cokernel vector of C{6i,dj){s). By (3.4), follow the relations 

%,{s)A{S,){s) = a(g'^-i.)3ff^S,,-(f ), 

Vi,j{qs) 

^For d = 2, the v(s), v{qs) in (3.5) here correspond to the gi,gi in [5] (9.8.17). 



'±=pk^-sq^)^'' 

\/n 



d-1 
k 
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where v,,,(s)'^-^=%(5)(e'^-i). 

We now consider the case of 5i, 5j with the relation Sj = q'^^^'^^Si. Choose the diagonal matrix 
M in (3.1) as 

5 = dia[l,g'^-\---,g(^-i)('^-i)], (3.9) 

i.e. Si = g^''"^)' for < i < AT — 1. As the transfer matrix t(s) (of even chain-size) is the same 
when changing the L-operator L(s) in (2.1) to L(— s), we shall use Sjj = Sij{—s) in (3.6) for 
the construction of Q/j-operator, (i.e. replacing a{s) by o(— s), and s by —s for the vector Vj j in 



-(d-l) 



1. T(rf-l) 



(3.5)), but keep Sij = Sij(s) in (3.8) for the Q^-operator. By setting {vij)o = Sj ^ {sq 2 ) 

(vjj)^ = (5j ^ {sq~) 2 for vectors in (3.5), Sjj, Sjj are defined to be zero except i—j = ±1 G Zj^, 
in which cases the in (3.6), (3.8) are given by^ 



v.,,(=v,,,(s)) = ||g('^-i)^-(sg^)^-), 
Vij(= v„-(s)) = (g-('^-l)^(sg^)^--'^| 

Here the vector \\s) and covector (s|| for s G C* are defined by 



d-l 



fc=0 



d-l 
k 



k=0 



d* 



The S-operator is now in the form 

/ So,i 
Si,o 



s = 





Si,2 











So, 



Af-1 



\ Sjv-1,0 



Sjv-i,Ar-2 





SAr-2,Ar-i 




and the same for the S-operator. By (3.7), one finds 



A{6j){-s)Sij{s) = a{-s)q 2 Sij{qs), 

2 Sij{s/q), 



D{di){-s)Sijis)=ai-q'^-^s)q 
then by the even chain-size, follows 

t{s)QR{s) = a{q''-h)'^QR{sq-^)+a{s)''QR{sq). 

Similarly, one has 

QL{s)t{s) = a{q'^-^sfQL{sq-^)+a{s)''QL{sq). 
In order to construct the Q-operator from Qr and Ql, it suffices to show ([3] (C28)): 

Ql{s')Qr{s) = Ql{s)Qr{s') for s, s' G C*. 



(3.10) 



(3.11) 
(3.12) 

(3.13) 



^Note that for d = 2, vectors in (3.10) here are the same as [29] (4.4) and (4.7). 
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To do this, as in [8, 29] we consider the product function of (3.10) , f{s', s\i,j; k, I) := Vjj(s')vfe,/(s), 
which is equal to zero except \j — i\ = \l — k\ = 1, in which case f{s', s\i,j; k, I) is expressed by 



r 



(t-i)(d-l)((d-l)-2r) 

d 

r=0 

Next look for an auxiliary function P{s', s\n) for n G Zjv such that 

fis, s'\i,j; k, I) = P{s', s\k - i)f(s', s\i, j- k, l)Pis',s\l - (3.14) 

by which the product Ql{s)Qr{s') differs by the boundary contribution when interchanging s by 
s', 

Ql{s)Qr{s) = P(s', s\ki - ii)QL{s')QR{s)P{s', s\kL+i - «l+i)~\ 

hence follows the relation (3.13) due to the periodicity of boundary condition. There are four cases, 
\j — i\ = \l — k\ = 1, in (3.14) to consider the function P, which automatically holds for j — i = I — k. 
The other two cases yield just one condition on the function P: ^p(gf = F(g(/5}nj ' where F{s\n) 
are defined by 

F(c,|n) = ^\(-+i)(<i-i)(^+'-)5^+r [<^-i] , (sG C*,nGZjv). 

r=0 ^ 

Since F{s\n) = F{s~^\n') for n + n' + 2 = (mod N), there are functions P{s', s\n) satisfying the 
relation (3.14). Hence follows (3.13). 
We now define 

Q{s) := Qr{s)Qr{so)-^ = Ql{sq)-^Ql{s), (3.15) 

where sq is a fixed value of s so that Qr{sq) and Ql{sq) are non-singular. By (3.11), (3.13), one 
finds [Q{s),Q{s')] = [t{s),Q{s')] = 0, and the TQ-relation 

t(s)Q(s) = a{q'^-h)^Q{sq-^) + a{s)^Q{sq). (3.16) 

As in the discussion of Q-operator of the eight- vertex model in [3, 19], the non-singular property 
of Q r{so) , Q l{so) for a general sq with generically arbitrary Tij,Tij in (3.6), (3.8) is expected, 
(unfortunately I know of no simple way to prove it). Nevertheless we shall assume the non-singular 

property of Qr{sci),Ql{so) for some sq, and define the Q-operator in (3.15), which is independent 
to the choice of Tij, Tij regardless of Qr, Qls dependence on them. Since vectors Vjj, Vjj in (3.10) 
satisfy the relation, g~('^~^)^^Vjj(s) = Vi+ij+i(s), Vjj(s)gf*^'^~^''^^ = Vi+ij+i(s), the relations 

will hold if the parameters Tij,Tij in (3.6), (3.8) are chosen with Tjjg'-'^"^^^^ = Tj+ij+i, q~^'^~^^^''Ti,j = 
Tj+ij+i, (e.g., by setting Tjj = yij{c), Tij = yi,j{c) for some fixed complex number c). This in turn 
yields Qr, Ql, hence the Q-operator, commute with q('^~^)^^ . When d — 1 is relatively prime to N, 
e.g. the cases d = 2, 3, — 1, one has [Q{s),q^''] = 0. 
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For j > 2, the jth fused L-operator for the XXZ spin chain is the matrix L^^\s) = 

(j'kli^))^^!^^^- C-'-auxihary and C^-quantum space constructed from (2.1). The C^-auxihary 

space is the space of completely symmetric {j — l)-tensors of with the basis e^"'^ = x^~^'~^y^ 
and the dual basis e^^* = {^~]r'^)x^~^~'^y'^ {k = 0, . . . ,j — 1), where x,y denote the standard 
basis of C^-auxiliary space in (2.1), x,y the dual basis, the symmetric (m + n)-tensor of 

defined by 



X y =x®...0x^y^...^y + others by permutations, 
\ n I ^ V ' ^ , " 

^ ' TO n 

and the same for x'^y'^. The C^-operator L^j^j(s) is defined in [30] with the expression 
^g(^) = Ji'h-f'^l'K ietU^) ^5aux l{qs) ®aux • • • ^aux l{q'-'s)\e^^) . 



The _7th fusion matrix is the trace of the monodromy matrix 

L 

I 



t(^)(s) = tr(^,((g)4^\s)), L'f\s) = L^^\s) at site £, (3.17) 



which form a family of commuting operators of (8) C^. Since L^'^\s) = Ylf=i a{q^ s)L{s) , we have 
r(2)(s) = l\tZiaiq's)h{s), hence the TQ-relation by (3.16): 

T(2)(s)Q(s) = h{s)^Q{q-h) + h{q-^sfQ{qs) (3.18) 

where the function h{s) is defined by 

d-l 

h{s){= ha{s)) = ^a{q's). 

1=1 

Using [30] (4.18), one finds the fusion relations among T^^^'s by setting T(°) = 0,T(^) = h{q-^s)^: 

rO)(s)r(2)(gJ-is) = h{q^-h)^T^^-^\s) + h{q^-hfT^^+^\s) (3.19) 

for j > I. By the induction argument, the TQ- and fusion relations, (3.18) , (3.19), in turn yield 
the r(-')Q-relation 

for j > 1. Then follows the boundary fusion relation: 

T{N+i) ^ j,{N-i) (g^) ^ 2h{q-h)^. (3.21) 



^Note that -L^ ;(s) here differs from \-i^]{s) in [30] by a factor, L'j^lis) = ^k^ll" " \-kl{s), contributed from 



the normahzed factor in the (corrected) formula (4.27) of [30]: T^^'{i) := -p-rw+j-a r — z — -ij > 1- Hence the 



LLk = d-2 ^ ' 



fusion matrix T'-^-'(s) in (3.17) of this paper and t^^^is) in [30] are related by T^^\s) = ^^T/_2 1 

lli=d-i "'^"''^'^ 



is). 
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Similar to the functional equation satisfied by the chiral Potts transfer matrix ([9] (4.40)) in 
CPM , the Q-operator of XXZ spin chain of higher spin encoding the root-of-unity property is 
expected to satisfy the Q-functional equation as in the cases of six-vertex and eight-vertex models 
[19, 20, 29, 31]: 

where the s-automorphism C is conjccturally defined by C{s) = —s, and Mq is some normalized 
matrix. By (3.20), the Q-functional equation (3.22) is equivalent to the A^th QQ-relation 

T^^\qv) = M^^Q{Cs)Q{s), (3.23) 

or either one of the following QQ-relations: 

rO)(gs) + T(^-^){qJ+^s) = M^^Q{Cqh)Q{s) (3.24) 

for < j < N. Note that from the structure of eigenvalues of the Q-operator, the TQ-relation 
(3.18) in turn yields the well-known Bethe equation of spin-^^ XXZ chain (see, e.g. [22, 30] and 
references therein): 

with the evaluation parameters of the s/2-loop-algebra symmetry (see [30] (4.32)) built into the 
QQ-relation (3.23). 



4 The Matrix Q{s) of XXZ Spin Chain of Higher Spin for a Gen- 
eral q 

We now study the XXZ spin chain of higher spin with an arbitrary value of q (not necessary a root 
of unity), which we will assume in this section. By imitating Baxter's method of producing the 
Q-operator of six- vertex model in [5] section 9.8., we construct a Q-operator of XXZ spin chain of 
spin {d > 2) for a general q. 

We search vectors as columns of Qr{s), whose image under t(s), by (3.11), is the sum of two 
proportional vectors with twisted variables. To do this, we perform the gauge transform of the 
monodromy matrix (2.4) by choosing two-by-two matrices Pi, . . . , and Pl+i = Pi of the form 

^^"(r^ 1) (l^^^^ + l)' ^L+i=n, 

such that the bottom-left entry of P^~^L^(s)P^+i has a non-trivial kernel vector in C*^ for 1 < ^ < L. 
Under the gauge transform by P^'s, the L-operator (2.1) at site i of the monodromy matrix is 
replaced by 

P-L(.)P,,,= f ] 
' ^ ^ V C(r,,r,+i)(s), D{n){s) ) 
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where A(r£+i), D(r^) and C(r^, r^+i) are the operators defined in (3.2). By Lemma 3.1, the condition 
of C{re,ri+i) with the non-trivial kernel is given by r^+i = r(,q'^~^~'^^ for some < A; < d — 1. Thus 

the general solution for r^-'s is : = 7.(jr^i+-- +^'?-i , where r is arbitrary and each takes the value 
d — \ — 2k for < < d — 1. The boundary condition rj^+i = ri is satisfied if Si + . . . + ei, = 0. 
Let V£(= V£(s)) be the kernel vector of C(rf,r£+i) in Lemma 3.1 with {mi)o = 1. The product- 
vector (8'^^]^v^(s), depending on r and Si, . . . , e^, is denoted by v(s; r, e)(= v(s; r, ei, . . . , El))- Using 
the first two relations in (3.4), one finds the vector v(s;r, e) under the transfer matrix (2.5) is 
decomposed as the sum 

t(s)v(s; r, e) = a{s)^\/{qs; r, e) -|- a{q'^~^ s)^y{q~^ s; r, e). 

In later discussion, we need the explicit form of v(s; r, e). For this purpose, we restrict the selection of 
£ only in the form = ai{d— 1) with = ±1, and write y{s\r,a){= y{s\r,ai, . . . jCFl)) = \/{s;r,s). 
Now the re is 

= rqi'i-m'^i+-+'^t-i) ^ with ct^ = ±1, cti + . . . + ctl = 0, (4.1) 
a condition consistent with even L. By (3.5), the product-vector y{s\r,a) is expressed by 

y{s\r,a) = <S> [ r^'q^'^'^~^^^i=^'''{-sq-^Y'^' ^ke] 4' ' (4-2) 
(.=1 \kc=0 * / 

and satisfies the relation 

t(s)y(s|r,cr) = a{s)^y{qs\r,a) + a{q'^''^ s)^y{q~'^ s\r,a). 

The operator Qr{s) is a (d — 1)^ by {d — 1)^ matrix whose columns are linear combinations of 
vectors y{s\r,a). Hence follows the relation (3.11). Similarly, using the covector-relation in (3.4), 
(3.5) of Lemma 3.1, one finds the product-covectors 

y{s\r,a) = (g) ( E r-^^^g-^^^^^-^) Eti -^(sg^)-^'=^efe,,, ) , (4.3) 
t=i \ki=i) J 

satisfying the relation, 

y(s|r, cr)t(s) = a{(f~^ sfy{(f^ s\r, a) + a(s)'^y(gs|r, cj). 

Hence the relation (3.12) holds for QL^operator constructed from the covectors y(s|r, cj)'s. 

To construct the Q-operator, we need the commutative relation (3.13) about Qr,Ql- To do 
this, we consider the scalar product y{s'\r,d)y{s\r^a) of two vectors in (4.2), (4.3), now evaluated 
as 

L / d-1 _ ^ _ ~ _ \ 

Jl ^(_I)fe« r^-ij g'=^(<^-i)(Ei=i^<-Ei=i<^<)(sg^)<^<'=^(sg^)^«'=« . (4.4) 



As the argument in [5] (9.8.23), we show the above product function is a symmetric function of 

y{s'\r, a)y{s\r, a) = y{s\r, a)y{s'\r, a). (4.5) 



s', s: 
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It is obvious that (4.5) is valid if ct^ = cr^ for all i. For a given a, all values of 5^ allowed by (4.1) are 
obtained by successive interchanges of pairs (cfj, aj+i). It suffices to show the symmetric property 
of the ratio function 

(y(s'|r, ...aj+i,aj, ...)y(s|r, a)) / {y{s'\r, ...aj,aj+i, ...)y{s\r, a)) . 

Since all terms but £ = j,j + 1 in (4.4) are symmetric in aj, cfj+i, the above ratio leaves only these 

terms in numerator and denominator. The direct calculation on cases for aj,aj^i,aj,aj+i = ±1 

reveals the ratio function is indeed a symmetric function of s' , s as the case d = 2 in [5]. Thus one 

obtains (4.5). As any column of Qr{s) is a linear combination of vectors y{s\r,a) and any row of 

Ql{s) is a linear combination of covectors y{s\r,a), (4.5) in turn yields (3.13). We now consider 

the set of vectors y(s|r, a) formed by all possible complex numbers r, and all possible (^^2) values 

of cr = (c7i, . . . ,c7l) in (4.1); similarly for the covectors y{s\r,a). As in the six-vertex model case 

([5] page 198), conjecturally there are values of s for which these vectors span the quantum vectors 
L , L , 

(8) C , (8) C respectively, so that Qr{s) and Ql{s) can be chosen non-singular. Let sq be a such 
value and we defined the Q-operator by (3.15). Then the matrix Q{s) satisfies the TQ-relation 
(3.16) . 

5 Concluding Remarks 

In this work, we first construct a Q-operator of XXZ chain of higher spin for the root-of-unity q, 
then another Q for the general by the similar method of producing the six-vertex Q-operators 
in [29, 5], respectively. These two constructions provide different (though related) Q-operators of 
XXZ chain of higher spin, with the same TQ-relation. The methods may help us to understand 
better the Q-operators of solvable lattice models in general. Usually Q-operators in root-of-unity 
cases are often of particular interest, as shown in the pure ice model. In the XXZ chain of higher 
spin at root-of-unity g, the Q-operator with the constraint of Q-functional relation (3.22) is the 
most important one as it encodes the symmetry of the system. Unfortunately at present I don't 
know any way to prove the complete set of functional relation using the Q-operators in this work, 
despite the fact that with a similar Q-operator, the functional relations in the spin-^ case have been 
justified in [29]. For higher spin cases, a better quantitative understanding about the Q- and fusion 
matrices should be essential and required for further development on the root-of-unity symmetry of 
the theory. A programme along this line is under consideration, and progress would be expected. 
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